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Abstract 

We present a concept for on-orbit radio detection of ultra-high energy cosmic rays (UHECRs) 
that has the potential to provide collection rates of ~100 events per year for energies above 
10^" eV. The synoptic wideband orbiting radio detector (SWORD) mission's high event statistics 
at these energies combined with the pointing capabilities of a space-borne antenna array could 
enable charged particle astronomy. The detector concept is based on ANITA's successful de- 
tection UHECRs where the geosynchrotron radio signal produced by the extended air shower is 
reflected off the Earth's surface and detected in flight. 
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1. Introduction 

The origin of ultra-high energy cosmic rays (UHECRs) remains a mystery decades after their 
discovery |lll|2l|3i|. Deflection due to cosmic magnetic fields impedes the reconstruction of source 
direction except for the highest energies where current observations are statistically limited. A 
key to identifying the sources of UHECRs is to increase the exposure at energies above 6 x 
10'"^ eV. We present a space-based mission concept using radio detection techniques with the 
potential to increase the exposure by factors of 40-70 over Auger observatory at super-GZK 
energies relevant to charged particle astronomy. 

Radio detection efforts for UHECRs date to the 1960's fT] and remained active through the 
1970's 15, 6, 7, 8, 9, The technique originally promised an alternative to particle 

counters on the ground with a wider field of view. However, the effort was abandoned due to 
technological limitations. The field has had a relatively recent rebirth fTT] with the advent of 
low power fast digitizers. Existing ground arrays are currently adding improved radio detection 



capabilities lll5Lll6ll . 

In a different radio detection effort, the ANITA experiment was developed as a balloon-borne 
antenna array designed to detect radio pulses in the 200-1200 MHz frequency range from ultra- 
high energy neutrinos interacting in the Antarctic ice iflTl . No neutrinos were detected but there 
was a serendipitous observation of UHECR pulses reflected off the ice-sheet iTs*]. 

The Synoptic Wideband Orbiting Radio Detector (SWORD) is a mission based on the ANITA 



detection of UHECRs iH]. SWORD detects the geosynchrotron radiation ([llQ of UHECR 
extended air showers (EAS) produced in the atmosphere, reflected off the Earth's surface, and 
detected from low Earth orbit (LEO) (see Figure [!]• Orbit altitudes (600-800 km) provide a 
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Figure 1 : Diagram of the SWORD detection concept. The cosmic ray comes in from outer space and interacts with the 
atmosphere to produce an extended air shower The geomagnetic field induces a charge separation to produce a forwai'd 
directed geosynchrotron pulse which illuminates a spot on the surface of the Earth. The reflected pulse propagates though 
the atmosphere and the ionosphere to be detected by an antenna in low Earth orbit. 

large collection area for UHECRs with energies > 10^° eV with detection rates of ~ 100 per year 
However, these altitudes also require that the cosmic ray induced radio pulse propagate through 
the ionosphere, which is a major technological challenge for the success of this mission and is 
addressed in detail in this study. 

In this paper we describe the concept and analysis of SWORD. Section 2 presents the signal 
simulations based on a geosynchrotron radiation parameterization including models for oblique 
reflections, surface roughness, noise backgrounds, and dispersion the radio impulse through the 
ionosphere. Section 3 describes the mission architecture which entails the choice of antenna and 
triggering of pulses dispersed through the ionosphere. Section 4 gives the expected performance 
of the mission based on detailed simulations. In Section 5 we give the conclusions of our results 
and the outlook for developing this mission. 



2. Signal Simulation 

In this section we provide a parametrized model of the signal production from geosynchrotron 
radiation with propagation effects due to reflection on the surface of the Earth, including surface 
roughness, and the ionosphere. The noise contributions of the sky, ground, and anthropogenic 
backgrounds are estimated. 

2.1. Geosynchrotron Radiation 

Geosynchrotron radiation is produced by the electrons and positrons in the extended air shower 
interacting with the geomagnetic field. The Lorentz force separates the charges causing an overall 
transverse current [19, 20|. The total number of electrons and positrons in the shower propagates 
in a meter-scale longitudinal thickness 'pancake' resulting in a radio pulse that is coherent at 
frequencies below 100 MHz with partial coherence up to <1 GHz \ At these frequencies, 
the radio impulse power is proportional to the square of the shower energy resulting in strongly 
detectable pulses over large distances at the highest energies. 

The parametrization of geosynchrotron radiation applied in this study is based on the results 



of a data-driven maximum-likelihood model applied to the ANITA UHECR events ||18|] . Despite 
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steady progress in radio emission calculations over the last decade 112 iL l22l |23L l24l |25L 12611 these 
models do not directly apply to the geometry of a detector in flight where observations are in the 
far-field, in contrast to ground arrays. Observations from the sky also involve much larger zenith 
angles than ground arrays usually observe. The parametrization applied to the ANITA data and 
for the SWORD simulations is given by 
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The reference amplitude used in this parameterization is determined from the ANITA data 
at 360:i°° /iV/m/MHz 111. The values of Ao for parameterizations similar to that presented 
here vary widely in the literature. For our energy scaling with respect to 10''' eV, corrected for 
the exponentially falling spectrum, Ao ranges from as low as 30 juV/m/MHz j,2 7.] to as high as 
390 juV/m/MHz i2^ with values in between lying closer to the smaller value |ll3,|21, 29ll . 



The Rref I {R + Rxmax) term accounts for the l/R propagation loss of a radiative electric field 
including the distance Rxmax from shower maximum to the ground and the distance R from the 
ground to the detector The reference distance Rref - 8 km corresponds to the distance from 
shower maximum to the ground for a shower with 0. - 60°. 

The beam width of the UHECR emission is based on the synchrotron radiation formula found 
in ll30ll with parameters fitted to the ANITA data ifisll . The beam pattern parametrization F{f, Osh) 
depends on the frequency /, the angle of observation with respect to shower axis 0,/,, and is 
normalized so that its peak value is unity. The formula given in 13 Oil for the intensity frequency 
spectrum for a differential soUd angle dQ. is 
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c is the speed of light in vacuum, y is the Lorentz factor, is the angle of observation (in this 
case 6si,), p is the radius of curvature of the particle, and K2/3i^) and Ki/i(^) are modified Bessel 
functions. The beam parameterization is normalized to the spectrum at = and given by 
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where n is an exponent that is fit to match the ANITA data. Inserting Equation|2]into Equation]!] 
reduces F to 
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where ^0 - cop/{3cy^) 

The radius of curvature for an electron or position under the influence of the geomagnetic field 
B ~ 50 fiT with y ~ 60 gives p » mcy/3/eB « 2 km (assuming a 7=60) 1 1 81I . At frequencies 
between 30 MHz and 1 GHz and at observation angles of order 6 ~ 1°, ^ ranges between 6 x 10 
and 10^5. If ^ « 1 then Ki/3(^) <sc K2/3{^) giving 
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With the values of 7 ~ 60 and B ~ 50 yuT, we have the parameterization 
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where n - 3.333 is the exponent that best fits the ANITA data. Figure |2] shows the beam pat- 
terns as a function of Osh for selected frequencies as well as the half power beam width of the 
parameterization as a function of frequency over the range 30-300 MHz. 

The spectral shape 5'(/) follows an exponential drop according to the ANITA observations 
in the range of 300-1000 MHz lii. The extension to lower frequencies is guided by modern 



simulations 11251 12611 that predict that the spectrum of inclined showers, which SWORD is most 
sensitive to, flattens out below 100 MHz. It is unclear whether this feature remains for far-field 
observations such as ANITA or SWORD but the spectral flattening assumed in the simulations 
is a conservative approach and the results would only improve if the signal strength continues to 
rise with decreasing frequency. 

The Bj^ term is proportional to the cross product of the shower axis with the geomagnetic field. 
The shower strength varies with the magnetic field intensity, which is referenced to Bref = 45 fiT. 



2.2. Oblique Reflections 

The remaining terms of the parameterization in Equation [1] account for reflection effects on 
the surface. The reflection is calculated using the scalar Kirchoff theory of diffraction. With this 
formalism we estimate the frequency and distance dependent coefficients for oblique reflection 
of radio signals off the surface of the Earth. The relevance of this approach will become evident 
in the treatment of surface roughness. 
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Figure 2: The geosynchrotron beam pattern for various frequencies. Shown on top are the F{f, 8si,) factors in Equation^ 
which gives the beam pattern as a function of angle with respect to the shower axis. The function is nonnalized to the 
peak emission. The plot on bottom shows the half-power beam width (HPBW) as a function of frequency. 



The geometry of the oblique reflection is shown in Figure |3] The source is located at point S 
and the receiver at point P. The origin is located at the specular reflection point and the position 
vectors of the source and payload are Ri and R2, respectively. The zenith angles for each position 
vector are 6^ given by the specular reflection condition in the plane of incidence. The position 
vector p-(p,(f)) of an area element on the surface is dA. The vector from this element to the 
source is r' = Ri - p and to the receiver is r = R - p. In the Kirchoff theory of diffraction, the 
signal at point P has contributions from reflections of the source at S over the whole surface. 

The region over which the path length difference is less than some fixed value will be bounded 
by an ellipse on a flat surface, with major radius and minor radius b, satisfying the condition 
a = bj cos 9^. Figure |3] shows an annular region with thickness dp. The problem as formulated 
so far is the same as determining the size of that region in a mirror around the specular point 
which contributes to the far-field reflection. In our case, however, the reflection coefficient is not 
that of a perfect conductor, so we must use the Fresnel coefficients for reflection from a dielectric 
surface. 

For this analysis we assume that p <s: Ri,R2 so that the polarization vector is approximately 
constant across the reflective region, and we can use a constant Fresnel reflection coefficient for 
each component of the polarization. Using the Kirchoff approximation |30], the received field 
strength for monochromatic waves at P due to a source field strength Esrc at point S is given 
by 

2^«Jw/fl« |r| |r'| 

where k - 2n/A is the wavenumber, the index p refers to the polarization vector component in the 
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Figure 3: Geometry for the reflections calculated using Equation |9] This is an adaptation of the standai'd Kircholf- 
Fresnel diffraction through and opening in an opaque plane. The illumination from source S has a projected ellipse on 
the ground which is reflected to a receiver at point P. The Kirchoff-Fresnel integral includes contributions not just from 
the specular reflection path, given by vectors Ri and R2, but from all paths on the projected area given by vectors r and 
r' . See text for details. 



transverse electric or transverse magnetic polarization, and the corresponding Fresnel coefficient 
'Fp(0) for a given zenith angle cos 6* = z ■ r' at dA. Let us ignore the effects of a curved Earth for 
now and assume that the reflection surface is planar 

As shown in Figure[3]we evaluate the integral over the projected ellipse of the source since this 
is the shape the UHECR radio illumination takes on the ground. For an ellipse with semi-minor 
radius b, given by + cos^ 0~ = b^, Equation[8]becomes 



^ /^-i/cosfi- nyjb--y^- cos- 8- ^ik\r\ ^ik\r'\ 

E'rcv^^.] dy \ ^ af^£L(w,0)5>(0)— — COS0 (9) 

J -hi CO!, 6, J- Vi2-y2 cos^flj |r| |r I 

Let us assume that b is small enough so that 6 « 6. over the region of integration and the 
half -power beam- width (HPB W) of the UHECR beam is approximately constant. The distances 
r ^ R\ and r' x R2 in the denominator of the integrand but not in the argument of the phasors. 
Equation|9]becomes 

k cosfl r-b/ cone, oyfc2-/cos2ft 

£L = £L(o.)W.)f^^ dy , (10) 

Zm K1K2 J-i/cosfl; J-^b^-y^cos^e, 

In this situation, the calculation reduces to an area integral weighted by the product of the two 
phasors e'*'""' and e'*^'"" '. In the case of propagation of a free spherical wave, the equivalent integral 
gives approximately the area of the Fresnel patch [31] 

AR2R\ 

Af - . (11) 

R1+R2 

We can estimate the degree to which the area integral approximates A/? by evaluating the integral 
in Equation [To] numerically. Figure |4] shows the numerical evaluation of the integral in Equa- 
tion[TOl scaled by the area of the first Fresnel zone projected on the ground Afl cos 9^, evaluated 
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at the Fresnel distance Rf - ■\/ARiR2/{R\ + ^2)- The 1 / cos 6- factor is due to the obHque pro- 
jection to the elHpse. For high zenith angles, the integral converges to Af/ cos6- after the first 
couple of Fresnel distances Rf. The integral oscillates more as the zenith angle is decreased to 
vertical where the source distance to the reflection point is smaller (see Figure |4}. The distance 
from the specular reflection point to the UHECR shower maximum, assuming an average shower 
maximum depth of 725 g/cm^, and to the satellite, assuming an altitude of 800 km, are shown 
in Figure |5] For zenith angles above 30°, which are of interest to SWORD, it is safe to assume 
that the magnitude of the integral is given by Af/ cos6~. The magnitude of Equation [TOl then 
becomes 

£L * EUco)p^, (12) 
Ri + Ri 

as expected from far-field specular reflection. The oscillations shown in Figure|4]are of the order 
of 10% after a few Fresnel distances. These oscillations are less pronounced for high bandwidth 
signals such as UHECR pulses compared to the monochromatic estimates shown. 




Figure 4: The area integral of Equation llOl evaluated for variou.s zenith angles. Equation llOl is an approximation of the 
Kirchoff-Fresnel surface integral of Equation |8]where the integrand consists of the phasors corresponding to the different 
path lengths light can travel. The result shows that only the first Fresnel zone, projected on the ground, contributes to the 
magnitude of the integral. The oscillations due to diffraction become more stable over a wider region of integration. 

The last piece to justifying these approximations is to compare the size of the UHECR radio 
illumination, projected on the surface of the Earth, to the size of the Fresnel zone. Figure |6]shows 
that ratio of the projected UHECR beam HPBW ellipse, for both the semi-major and semi-minor 
axes, to the Fresnel distance Rf. The size of the ellipse is in all cases greater than Rf. For 
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projected HPBW to Rp ratios greater than 3, where the area integral converges, the assumption 
that Esrc is constant in Equation|9]is well justified. Greater care may be needed for ratios lower 
than three but, as will be shown, there are not a lot of events in the SWORD field of view for the 
low zenith angles where this is the case. 




Figure 5: Distance to tlie speculai' reflection point as a function of zenitli angle for the cosmic ray shower maximum and 
for a satelUte at 800 km orbit altitude. 

The size of Fresnel zones Rf, for wavelengths between 1 m and 10 m has a maximum of 
~2 km for the geometries of interest. Figure Q shows the deviation from flatness due to the 
Earth's curvature for a patch of size Rf scaled with wavelength. The deviations from flatness 
are well below one wavelength for the range of interest to SWORD. These corrections may be 
of interest to an energy reconstruction analysis but for the purposes of simulating the experiment 
they are expected to be negligible. 

2.3. Surface Roughness 

Understanding how the reflection is affected by surface characteristics is crucial to the de- 
tectability and energy reconstruction uncertainties of the showers observed by SWORD. Surface 
roughness is incorporated into the reflection under the Kirchoff approximation of scalar fields, 
treated above, via the Rayleigh criterion for the diffraction effect of the roughness. The rough- 
ness model implemented in this study follows the self-affine fractal models of the surface of the 
Earth used in radar reflection analysis lls^l . 

Many, if not most, topographic surfaces display fractal behavior which is not strictly self- 
similar, but which has a self-affinity such that the surface roughness scales with the distance over 
which it is evaluated |32]. The root-mean-square (RMS) roughness cr/,(L) - V<~Cz--'iP^i for 
measured heights z, scales as a power of the distance L over which it is calculated. The depen- 
dence can be parameterized in terms of a reference scale Lq and the corresponding roughness 
value at that scale is ^ 

crh(L)^crh(Lo)l^] , (13) 




Zenith Angle (degrees) 

Figure 6: Ratio of the projected UHECR beam HPBW (see Figure |2) ellipse semi-major and semi-minor axes to the 
Fresnel distance Rf = ^jAR[R2/(Rl + Ri) calculated from the distances to the specular reflection point from shower 
maximum and from a receiver in orbit at 800 km altitude (see Figure [5). The horizontal dashed line marks a ratio of 
three, above which value the Kirchoff-Fresnel integral is stable (see Figure^}. 

where H is known as the Hurst parameter, Q < H < 1 . Typical values for H are of the order 0.5 
for many topographic surfaces. 

For the SWORD simulation we treat three types of surfaces: land, ice, and ocean. For land 
we assume a typical Hurst parameter of 0.5. Although the cr/,(l m) can vary greatly we assume 
a typical value of 20 cm. As we will see, it is not expected that many SWORD events come 
from reflections on land due to the large amounts of anthropogenic noise present. However, 
for analysis of the data these parameters will have to be better estimated. For the ice we have 
estimated cr/,(120 m)=0.05 m and a Hurst parameter of 0.65 from in-situ measurements at Taylor 
Dome, Antarctica |33]. An analysis, found in [33], of radarsat data |34] confirms that H ^ 0.65 
to length scales up to a kilometer For the ocean, analysis of significant wave height from Jason- 1 
data 1 35] gives cr/,(2 km)* 2.7 m for the southern ocean. The SWORD simulation uses a 30°x30° 
grid of cr/,(2 km) values ranging from 0.5-2.7 m and assumes a Hurst parameter of 0.5. 

For surface roughness with RMS deviation cr/,(p) that is less than a wavelength, the scattered 
field strength upon reflection form a roughened dielectric surface is given by ll32ll as 

E rough = E smooth exp [-2k^(rl(p) cos^ 0;) . (14) 

In a Huygens-Fresnel application, such as the one in the previous subsection, where each sur- 
face area element is treated as re-emitting the incident radiation, we can assign this scattering 
amplitude to each area element, and thus include it explicitly in the integral of Equation|8] Note 
here that cr/,(p) is assumed to have explicit radial dependence, thus each annulus in the Kirchoff 
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Figure 7: The deviation from flatness for a distance Fresnel distance Rf and 3Rf as a function of zenith angle given 
in wavelengths. The blue solid line is given by REanhi^ ~ Rf cos(9-)/REanh) and the red dashed line by RiartlA^ ~ 
3Rf cos{9-)/Rgar,l,). Figure |4] shows that above 3Rf the contributions to the reflection are of order 10%. This figure 
shows that the deviations from flatness are well below one wavelength for the range of interest to SWORD. These 
coiTections may be of interest to an energy reconstruction analysis but for the purposes of simulating the experiment they 
are neghgible. 



integral will have an assumed common RMS height variation, which changes with the radii of 
the annulus. 

Including the surface roughness to the Kirchoff integral in Equation[TO]gives the expression 

eL. = Kci(-^)rp{e,)——^ I dy \ dxe H>' v^^Vi'i (is) 

Zm K\K2 J-i/cosfl; J- sjh--y- cos- 0, 

The surface roughness is dominated by the Fresnel distance scale Rf- This means that the 
effective roughness contribution results in a roughness term 

G{f, R, 0,) ~ exp [-Ik^crliRp) cos^ 0,] (16) 

Examples of the attenuation due to surface roughness are given in Figure [8] The attenuation 
due to ocean surface roughness is of order a few dB for frequencies below 100 MHz and severely 
affects the signal strength at higher frequencies. For ice the attenuation is not significant in any 
of the frequencies of interest to SWORD. 

2.4. Dispersion of VHF Pulses Through the Ionosphere 

The ionospheric plasma has a well-known dispersive effect on VHF pulses with negligible 
attenuation. Extensive studies of dispersion were performed with the FORTE satellite 11361 137 , 



42l |44|1 . on which the simulations used for SWORD are based. When a VHF pulse propagates 



through the ionosphere, it is dispersed from several nanosecond duration up to 10-100 jUS severely 
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Figure 8: The roughness factor from Equation II 6 I for Antarctic ice and ocean. The ice roughness results in very small 
attenuation factors for the VHF frequency band. For the ocean, the attenuations factor is within several dB for frequencies 
< 100 MHz and becomes worsens significantly at higher frequencies. 



affecting the signal to noise ratio. This requires that SWORD dedisperse the ionospheric effects 
in real-time to restore the temporal coherence of the pulse and enable efficient triggering. 

There are several issues involved in developing an efficient trigger for space-based UHECR ra- 
dio detection. The time and direction of incidence of the radio pulse on the detector are unknown 
requiring the SWORD trigger to be continuously scanning for transients. The exact state of the 
ionosphere is not known a priori to the degree needed for efficient triggering, which requires the 
dedispersion algorithm to scan over a range of ionospheric transfer function parameters. The 
scanning range can be constrained by knowledge and predictions of the state of the ionosphere. 
The optimal trigger requires a matched filter implementation but, because the search space is 
large and spacecraft power is limited, we need to understand what is the minimum number of 
assumptions and resources required to perform the dedispersion efficiently. 

The implementation of ionospheric dispersion follows Issll and ll39[l and is described in the 
following. The index of refraction of the ionosphere is given by 

n\f) = (17) 



l_l(Mf]l^+ I UJWff sin' 13 (f ;f)?.i„2p 

^ 2 i-(fp/f) - V4 J > Sin P 



where = e^Neli^Ti^eome) is the plasma frequency, fc - eB/m,, is the cyclotron frequency, A^^ 
is the electron density, e is the absolute value of the electron charge, is the electron mass, eo 
is the permittivity of free space, B is the magnitude of the local magnetic field, and /3 is the angle 
between the geomagnetic field vector and the direction of propagation. The + in Equation [17] are 
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for the ordinary and extraordinary modes of propagation The frequency dependent phase shift 
due to ionospheric propagation is given by 



c Jt 



r)dr 



(18) 



where c is the speed of light, T is the location of the transmitter, R is the location of the receiver, 
and the integral is taken over the trajectory between T and R. 

Typical values are /c ~10 MHz and fp » 5 MHz in periods of low solar activity and fp x 
20 MHz for high solar activity. For the 30-300 MHz band we can perform a Taylor expansion of 
the index of refraction to give 



]^.f^ S fpfc 
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(19) 



where i = ±1 with a positive value for wave polarized opposite to the cyclotron rotation of 
electrons (ordinary modes) and negative for waves polarized with the cyclotron rotation (extra- 
ordinary modes). 

For a receiver at distance D away from the transmitter, the group delay is given by 



with the coefficients given by 
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An^ceonie Jt " In'ceoinl 

The total electron content (TEC) is the vertical column density of electrons between a trans 
mitter on the ground and a receiver a distance h directly above it given by 



TEC 



Jo 



(z)dz. 



(24) 



The TEC is given in units of TECU = 10'^ electi-ons/m 2. We define the slanted TEC (S TEC) 
for the case where the receiver is at zenith angle 0, with respect to the transmitter on the ground 
at a distance D. 
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For the practical implementation in the SWORD simulations we use the approximations used 
by Myre [39] giving S TEC = TEC/ cos 0„ 



C2 = aS TEC. 
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(26) 



with a = eV(87rVeom,) = 1.3445 x lO' Hz/TECU, 
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where /; is the receiver altitude and T is the thickness of the ionosphere for a slab model. 

The TEC values can be obtained on a daily basis from GPS derived data |40] and the magnetic 
field amplitudes and directions can be obtained from the international geomagnetic reference 
field (IGRF) model l4l[l . Figure |9] (top) shows the ionospheric distribution of GPS derived TEC 
values. The dispersion effects are expected to be best behaved near the poles and worst near 
the equator where the TEC has large values and fluctuations, especially when exposed to the 
sun. The stability of the ionosphere is dependent on solar activity meaning that for best results 
SWORD would need to be operate during a period away from solar maximum. 

As will be shown in Section 3, the real-time de-dispersion trigger is highly sensitive to the 
TEC. To save on operations count and spacecraft power it is important to constrain the TEC 
search space as much as possible. With the conservative assumption that the data can only be 
uploaded once a day, we estimated the prediction error between the ionospheric TECU over 
the SWORD sun-synchronous orbit by comparing the TEC values from one day to the next. 
Figure [10] shows the difference in TEC between two consecutive days as a function of orbit 
latitude. The largest deviations are near the equator, as expected. The overall RMS of the dif- 
ference is ~2 TECU as shown in the histogram in the bottom of Figure [10] A study with the 
FORTE satellite, comparing the reconstructed TEC using the observed dispersion of VHF wide- 
band pulses, compared to MAGIC, which is a GPS derived ionospheric solution, shows a scatter 
of a!2.8 TECU |42]. Altogether, it is expected we can bound the ionospheric TEC in the field of 
view of SWORD by about 3-4 TECU. 

2.5. Noise Backgrounds 

The dominant source of noise in the VHF band is the galactic noise background with an av- 
erage value of ~ 10"* K at 30 MHz decreasing as a power law with index -2.8 to ~ 10 K at 
300 MHz. To estimate the noise on SWORD we use the galactic background maps provided by 



1431. The SWORD antenna, which will be described later, has a 7 dBi gain. The noise maps of 
1 43]] are therefore averaged in 30° x 30° bins on the sky. This averaged galactic noise background 
can vary by about a factor of 8 depending on whether the galactic center is in the field of view. 

The SWORD antennas are pointed towards the Earth's horizon. This means that a portion 
of the galactic noise contribution is reflected off the surface of the Earth prior to exciting the 
antenna. For an antenna with gain g(f, 0, (p), the effective antenna temperature is given by 

Ldng{f,e,(/>)T,,,rceif,e,(/>) 

TaH) = J (29) 

For the SWORD geometry, the noise sources are the direct observation of the sky noise temper- 
ature Tsky, which consists of about 75% of the antenna beam pattern, and the observation of the 
ground, which consists of the remaining 25% of the beam pattern. The ground thermal noise is 
a combination of the reflected sky noise radiation, which is weighted by the Fresnel reflection 
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Time (s) 

Figure 9: (Top) distribution of the ionospliere at one point in time. The distribution of total electron content (TEC) 
given in TECU (lO'^ electrons/m^) is typically worse at equatorial latitudes during the day time. An 800 km altitude, 
80 degree inclination, sun synchronous orbit is superimposed on a continental outline. The small circle on the top right 
shows the size of the intersection of the ionosphere's peak density, at ~400 km, with the field of view of the payload. 
The larger circle shows the contour of the payload's horizon field of view. (Bottom) the TEC distribution directly below 
the payload as it orbits the Earth. 



coefficient T^Tsky and thermal emission from the surface (1 - T^W surface- The surface noise 
temperature is approximated at 240° for ice and 270° for land. The ocean is assumed to be 
completely reflective at VHF frequencies. 

The anthropogenic noise background is also modeled in the SWORD simulation. The FORTE 
satellite made extensive studies of the radio noise backgrounds on Earth |45, 4^. For the 



SWORD simulation, we summarize the effects of anthropogenic noise with a noise factor NF, 
which is a proportionality factor to the antenna noise temperature due to the sky noise and in- 
versely proportional to frequency 

/39 MHz\ 

Tan,h,o{D = TUfWF - 1) ( ^ 1 (30) 



The noise factor NF is strongest over North America, Europe, and East Asia, peaking at NF ~ 
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Figure 10: Difference in TEC (ATEC)from one day to the next for a sun-synclironous orbit. Shown on top is the ATEC 
directly below the payload for each of 16 orbits completed in one day as a function of orbit latitude. The histogram in 
the bottom is the distribution of the differences plotted above. 



40, and weakest over open ocean, where there is practically no contribution. 

The final source of thermal noise is the system temperature Tgys of the instrument itself. A 
typical value of Tgys = 90 K is used based on current amplifier technology. 

3. Mission Architecture 

The SWORD payload consists of an array of dual-polarized log-periodic dipole array anten- 
nas. Each antenna polarization has its own signal chain that is fed to a trigger system and data 
recording system. The payload can fit in a Pegasus Fairing envelope with 1.94 meter height and 
diameter that tapers from 1.16 - 0.72 meters. A version of the SWORD payload with 16 anten- 
nas is shown in Figure [TT] The orbit of choice is an 80° inclination, sun-synchronous orbit, with 
altitude ranging from 600-800 km. The high inclination is driven by observation of the poles 
where the ice provides a relatively smooth reflector with generally radio quiet locations. Placing 
the orbit at the terminator provides enough exposure to the sun to power the spacecraft while 
keeping away from the most excited portions of the ionosphere. 

The spacecraft will be equipped with a GPS receiver as well as attitude sensors for tilt and 
heading. The SWORD system diagram is shown in Figure [12] An additional processor will be 
used to discriminate against lightning and anthropogenic backgrounds. For signals digitized at 
600 Msa/s with 3-bit resolution and a record length of <100 fj.s in a payload consisting between 
6-16 dual-polarized antennas, the size of a raw event ranges between 2.2-5.7 Mbits. Based on 
satellite maps of lightning strike rates ghtning is expected to trigger the payload at average 

rates of ~1 Hz. Keeping the telemetry rates down to several Gbit/day will require efficient fil- 
tering of background events. Fortunately, it is well known that the UHECR signature is directly 
correlated to the direction of the geomagnetic field. Along with wide-band spectral and event 
isolation requirements, it is expected that the false positives rate will be low. 
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Figure 11: A 16 antenna version of the SWORD spacecraft. Two rings of 8 antennas each deploy from a Pegasus 
Fairing envelope. The cylindrically symmetric configuration provides synoptic observation of the Earth's surface, which 
provides a full-sky survey of UHECRs as well as an added handle for identifying radio frequency interference. 



3.1. Antenna 

The spatial constraints of the Pegasus envelope and VHF wavelengths, ranging from 1-10 
meters, requires that the SWORD antennas be deployable. The wire structure of the log-periodic 
dipole antenna naturally lends itself to this application. Deployable LPDA antennas have been 
flown on other missions, most notably FORTE. 

The design shown in Figure[T3]provides a radiation pattern with gain ranging form 7-5 dBi in 
the frequencies of interest. The HPBW of the beam pattern ranges from 100°-120° providing a 
wide field of view for UHECR detection. 



3.2. Trigger 

The ti-igger system uses the ANITA- 1 HIItI and ANITA-2 lEgl threshold and coincidence 
approach. The signal is split into bands to protect the system against RF interference. This way, 
bands can be masked if they are being saturated. Each band triggers on a low threshold at MHz 
rates with temporal coincidence and geometric requirements to bring the rate down to a 5-10 Hz 
payload trigger rate. In the case of SWORD, the individual bands have to be digitized in order to 
dedisperse the eff'ect of the ionosphere in real-time. 

The current trigger model uses 7 bands for each polarization with bandwidth distributed ac- 
cording to the the inverse frequency phase delay. The band limits are given by the recurring 
relation = /iA where A = (2,OQMHzl2>OMHz)^^\ 
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Figure 12: The SWORD system diagram. Each antenna polarization has its own signal chain consisting of passband 
filters and low noise amplifiers. The signals are fed into the central instrument with tunable notch filters for mitigation of 
caiTier wave interference. The second stage filter prevents any additional out of band noise picked up by the signal chain 
from being ahased. After a second stage of amplification the antenna signals are split into a trigger unit and a digital 
recording unit. 



/z.,„, (MHz) 


(MHz) 


Bandwidth (MHz) 


30 


42 


12 


42 


58 


16 


58 


80 


22 


80 


112 


32 


112 


155 


43 


155 


216 


61 


216 


300 


84 



Table 1: The SWORD trigger bands distributed according to inverse frequency phase delay. See text for details. 



The trigger requires that at least 4 of the 14 total bands (7 for each polarization) exceed a 
threshold of 1 .6cr where cr is the voltage RMS of the noise. In addition, the full band signal has 
to exceed 4.2cr. The threshold crossing accidentals rate is given by the cumulative distribution 
function of the Rician distribution for noise p = exp[-v^^^^^^y(2cr^)]. For a digitized signal, 
accidentals rate is given by the sampling rate multiplied by the probability of a threshold crossing. 
The threshold crossing accidentals trigger rate is estimated to be ~ 170 MHz for each sub- 
band while the full-band accidentals rate is r/i ^ 88 kHz. The probability that 4 out of 14 bands 
trigger is approximated by P ^'"^ C4(rshT)'^{l - rr)'" |48] where t is the trigger window. A 
T = 3 ns window gives P ^ 0.05. Combined with the probability that the full band is above 4.2cr 
(pfb - 2.6 X 10""^) gives a single antenna accidental probability of 1.3 x 10"^. At 600 Msa/s the 
total accidentals rate per antenna is » 1.1 Hz. For the full set of A^^ antennas, without geometric 
coincidence requirements, the total payload accidentals rate is given by A^^ x 1.1 Hz. 

Prior to applying the threshold and coincidence trigger, the signals need to be dedispersed. The 
algorithm applies coherent dedispersion via cross-correlation with a template. The dispersion 
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Figure 13: The SWORD 30-300 MHz log-periodic array (LPDA) antenna. The boom length is 4.5 meters long with 
longest dipole span of 3.7 meters. The wire structure of the LPDA lends itself to a deployable structure. 



phase, is given by 




(31) 



where C2 and C3 are the coefficients for the group delay (t = -d(p/da>) in Equation |20] The 
digitized voltage data for a given band is cross-correlated with the dispersion phase, with the 
opposite sign, to dedisperse the signal. The birefringence of the ionosphere acts on circularly 
polarized modes where the SWORD antenna observes linearly polarized modes. The template 
for each polarization is therefore derived starting from the assumption that the signal has po- 
larization angle 6p along Bj^ as defined in Equation [1] Let H and V define the horizontal 
and vertical polarization vector directions. The Fourier coefficient of the template is given by 
W = e'"'° (cos OpYl + sin OpY). The transformations between right/left circular and horizon- 
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tal/vertical polarizations are 

R = (H+ iV)l V2 



(32) 



L^(H- iV)l V2 (33) 
H^iL + R)/^/2 (34) 
V ^i{L-R)l ^ (35) 
The right and left circularly polarized modes of the pulse prior to dispersion are given by 

Wr = — (36) 

V2 

Wl = — (37) 

V2 

After dispersion through the ionosphere, Wr and Wi will carry additional dispersion phase. 

Wr = (38) 

V2 

Wl = p (39) 

V2 

Transforming back to H and V results in 



2 



2i 



(40) 



(41) 



An example of the birefringence of the ionosphere, shown in Figure[T4] is due to the combination 
of left and right dispersion modes. 

The dedispersion algorithm proceeds to cross-correlate the data with the template 

Uvpoi = Vvpot ® Wv (42) 

UhpoI = Vhp„i ® Wh (43) 

The parameters to which this algorithm is sensitive to are the TEC precision with which the 
dispersion can be modeled, the amplitude and direction of the geomagnetic field, and additional 
dispersion due to higher order terms in the expansion of the dispersion phase. The results of the 
trigger sensitivity to SWORD for these parameters is shown in Section 4. 
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Figure 14: Simulated cosmic ray event with 2 X 10"" eV propagated through an ionosphere with STECk9 TECU. The 
splitting at lower frequencies is due to the geomagnetic field. Although this pulse is strong enough to be observed above 
the noise background in this spectrogram, pulses at lower energies are buried in the noise. It is only after coherent 
dedispersion that the power at each band will be temporally coherent resulting in a pulse that can be distinguished from 
the noise background. 



3.3. On-board Data Filtering 



Limited telemetry rates require that SWORD have an efficient way to analyze data on board. 
Anthropogenic backgrounds are a major concern as they are likely to trigger the payload at a high 
rate. A spatio-temporal cluster rejection algorithm will need to be run on board similar to those 



successfully employed in the analyses of both ANITA flights Ill7ll48ll49ll . Strong anthropogenic 
backgrounds cover ~30% of the surface of the Earth. As will be shown in the next section, 
although the sensitivity of lO'^-lO^" eV is likely to be severely affected, events with energies 
> 10^" eV, which are most relevant for charged particle astronomy, are bright enough so that 
detectability will not be as strongly affected. 

In the remaining 70% of the surface of the Earth, mostly covered by ocean, lightning is a 
major concern. World lightning frequency maps (4T\ show lightning occurrence rates of 0.1-0.4 
events per km^ per yr over the souther ocean. For the SWORD payload, at an altitude of 800 km, 
with a visible area of 22 million km^, this translates to a rate of ~0.5 Hz. Although this seems 
like a low rate, assuming 3-bit digitizers running at 600 Msa/s with a record length of 100 fis (to 
accommodate for ionospheric dispersion) produces a data rate of 16 Gbits per antenna per day. 
With 12 antennas this is 190 Gbit/day. Data compression prior to transmission will be of aid but 
to keep telemetry rates reasonable the lightning false positives rate in the on-board analysis will 
have to be <10%. 

The main discriminator between lightning and geosynchrotron pulses is the requirement that 
the polarization of the UHECR be perpendicular to the geomagnetic field. Lightning is strong 
enough to pass the trigger thresholds on intensity alone but on-board analysis will be able to 
determine the polarization of the signal. The lightning pulse is also of longer duration than the 
UHECR and it is expected that better discriminators can be developed with further applications 
of the measurements from FORTE [sd] and the upcoming UHECR events from the third flight 
of ANITA 1,18.1. 
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4. Expected Performance 

The expected number of events accumulated by SWORD with a 3 year mission is shown in 
Figure [15] At energies beyond 3 x 10'^ the spectrum is extrapolated by either an exponential 
fall-off or a broken power law. The improvement factor over Auger, assuming a 2018 launch and 
3 year flight, ranges between factors of 40 to 70 at extremely high energies (see Figure [T6]l. The 
SWORD exposure for 10^*^ eV, compared to Auger and JEM-EUSO values from 12| is shown in 
Figure [TT] 




Figure 15: The number of events expected with SWORD assuming a 2018 launch and 3 years of operation. The top 
histogram shows the differential distribution for SWORD for both a power law and exponential extrapolation of the Auger 
spectrum. The event count spectrum for Auger, at the time of completion of the SWORD, is shown for comparison. The 
bottom histogram shows the expected cumulative event counts. 

The distribution of one year's worth of events are shown in Figure[T8]for energies > 3 x 10'^ eV 
and in Figure[T9]for energies > 10^*' eV. The absence of events on the continents is primarily due 
to anthropogenic noise. The concentration of events on the Antarctic ice cap with respect to the 
ocean is due to the smoother reflection surface of the continent. 

The distribution of zenith angles for SWORD triggered events is shown in Figure |20] At 
higher zenith angles, the SWORD payload has a larger collection area. However, the signal is 
attenuated at high zenith angles due to the obliquity factor in Equation[T] At lower zenith angles 
the conection area decreases resulting in a corresponding decrease in the number of events. 

Although an antenna array can reconstruct the direction of an incident plane wave with rel- 
atively high accuracy, as was demonstrated with ANITA [48J, the ability to reconstruct the di- 
rection of the UHECR event is limited by the unknown angle of observation with respect to the 
shower axis. The distribution of angles with respect to the shower axis for triggered events is 
shown if Figure|2T| The mean angle is 2.8° with a RMS of 1 .5°. It is expected that the interfero- 
metric direction reconstruction of the data will be significantly better than this so this distribution 
remains the dominant source of error for directional reconstruction. It is possible that radio spec- 
tral analysis could yield improved event direction reconstruction with application of techniques 

21 



Ratio of the SWORD and Auger cumulative distribution functions 
for power iaw and exponential extrapolations of the spectrum. 



100 



60 

S 



— Power law 

- - Exponential 














































































19.0 19^2 19.4 19.6 19.8 20^0 

loglO(Energy / eV) 



Figure 16: The improvement factor for the cumulative distributions shown in Figure [13] expected from SWORD over 
Auger assuming a 2018 launch and 3 years of operation. At super-GZK energies the improvement factor ranges from 40 
to 70. 
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Figure 17: The SWORD exposure compared to the Auger and JEM-EUSO values from [2]. The SWORD values shown 
here is for events with energy > 10^" eV for a 12 antenna version of SWORD. Unlike ground arrays, the exposure of 
SWORD is strongly energy dependent and rises steeply with increasing energy. 



similar to 115 111 . However, estimating the angle of observation with respect to shower axis using 
this technique does require that the attenuation due to surface roughness, which has a similar 
effect, be well understood. This direction is worth further exploration but we take the conserva- 
tive approach here and not include this in our estimates. The technique in jsill is currently being 
explored in simulations for application to the ANITA cosmic ray data expected at the beginning 
of 2014. 



The energy reconstruction is obtained by dividing the recorded voltage v/ at each frequency / 
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Global distribution for 1 year of events with E>3 x 




Figure 18: Global distribution for one year's worth of events with energy greater than 3 X 10 eV. The events are 
concentrated on the ice caps primarily due to the low surface roughness factors of ice (see Figure |8). The oceans are 
highly reflective but surface roughness is significantly larger, particularly at the higher end of the SWORD spectram 
(> 100 MHz). The low density of events over land is primarily due to increased anthropogenic backgrounds. 




Figure 19: Global distribution for one year's worth of events with energy greater than 10^" eV. Due to the large signal 
strength provided at these energies, the distribution of triggered events on the globe is not significantly affected by surface 
roughness or anthropogenic backgrounds. 

by 

Ef = (10'' eW)-^is(f)^ cos e,n9„f)] (44) 

where the direction reconstruction provides 9, and the geomagnetic field. The Fresnel coefficient 
Tifiz, f) can be estimated based on the surface of reflection (ice or ocean for most events). This 
estimate does not include surface roughness efifects or attenuation due to observing the UHECR 
signal off" the shower axis. 
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Figure 20: Distribution of zenitli angles observed by SWORD. The distribution peaks at around 70 degrees wliere the 
exposure is high but the signal is not attenuated by the obliquity factor. The decrease in the number of events at lower 
zenith angles is due to the reduced collection area corresponding to those directions. 
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Figure 2 1 : Distribution of observation angles with respect to the shower axis for the expected ensemble of triggered 
SWORD events. The mean value is 2.8 degrees with a root mean square of 1.5 degrees. The energy dependence of 
this angle tends to be wide for higher energies, where there is enough power the trigger away from the shower axis, and 
narrower for low energies, where the observation needs to be near the shower axis for the power to be above the noise. 



The total energy is estimated by a weighted average across the bands according to 

Y^fEjivfltiff 



-■measured 



Zf(vf/nf)^ 



(45) 



where is the noise RMS voltage. 

The distribution of reconstructed vs. true energy is shown in Figure|22] The increasing asym- 
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metry at higher energies is due to the fact that the UHECR radiation strength increases allowing 
the payload to trigger at larger angles off shower axis. In other words, the distribution shown in 
Figure|2T]is wider at higher energies and narrower at lower energies. Since we are not assuming 
that the angle of observation with respect to the shower axis can be known, it becomes a domi- 
nant source of uncertainty in the energy reconstruction as well as the direction reconstruction of 
events. The statistics of energy reconstruction are shown in Figure |23] The bias and errors are 
stable and well understood making a statistical energy reconstruction possible. 
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Figure 22: Distributions of the ratio of measured energy to true energy for lO", 3 X lO", and K)-" eV UHECRs. The 
width of the distribution is due to the combined effect of noise backgrounds and the distribution of angle of observation 
with respect to the shower axis, which is not accounted for in the reconstruction. 

The results shown so far do not take into account the dispersion effects of the ionosphere. In 
the VHP band the ionosphere does not significantly attenuate radio signals and the dispersion 
effects are understood well enough so that they can be removed in analysis. The greatest impact 
is on triggering efficiency. The ability of the SWORD mission to match the results shown in 
this paper depends on the real-time de-dispersion capabilities of the payload, which are limited 
by operations count and spacecraft power consumption. A full time-domain simulation of the 
dedispersion algorithm, using matched filtering for the ionospheric TEC and the geomagnetic 
field parameters was developed to guide the hardware requirements. 

The results of the simulations showing the trigger inefficiencies given the precision with which 
the ionospheric TEC is sampled is shown in Figure|24] To retain 50% efficiency at 5 x lO'*^ eV the 
dedispersion algorithm is required to sample templates in 0.05 TECU steps where the ionosphere 
can be constrained to within a few TECU (see Figure [TOli. 

A similar analysis to determine the sensitivity of dedispersion to uncertainties in the geomag- 
netic field amplitude and directions is shown in Figure |25] Knowledge of the geomagnetic field 
needs to be within 400 nT for 50% efficiency. The current IGRF models ll4l1l have errors of 10 nT 
making it more than adequate for dispersion template modeling. The stability of the ionosphere 
is critical for proper functioning of SWORD and periods of high solar activity, which can cause 
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Figure 23: Statistics of the ratio of measured energy to true energy. The mean values are error bars are derived from 
distributions such as those shown in Figure l22l 




Figure 24: Efficiency of real-time dedispersion trigger, relative to the trigger described in section 3, assuming cross- 
coiTelation template parameter errors in the TEC parameter. The SWORD trigger produces waveform templates based 
on a value for TEC and the geomagnetic field value. As the TEC value deviates from the true TEC, the efficiency drops. 



large fluctuations in TEC and magnetic field values, should be avoided. 

The duty cycle of the real-time dedispersion trigger is expected to be afi'ected by periods of 
high TEC. Figure |9] shows the TEC variation with orbit. The TEC rises to very high levels for 
10-25% of the orbit time. The SWORD payload would continue to operate in these intervals but 

26 



100 



80 



.1 60 
u 

LU 

I 40 

(U 
CL 

20 



10° 10^ 10^ 10^ lO" 10^ 

B-field Uncertainty (nT) 

Figure 25: Efficiency of real-time dedispersion trigger assuming deconvolution template parameter errors in the geomag- 
netic field parameter For this test the TEC is assumed true which the amplitude and direction of the geomagnetic field 
value deviates from the trae value. The cuiTent IGRF |41] model has geomagnetic field errors of 10 nT where the trigger 
efficiency is high. 

the trigger efficiency is expected to be affected. 



5. Outlook and Conclusions 



We have provided the scientific motivation, modefing, and expected performance for an or- 
biting UHECR detector using radio techniques that can provide the sensitivity required to open 
a new window to charged particle astronomy. Using the Earth's surface as a reflector achieves 
high exposure to energies > 10^" eV where cosmic rays are rigid enough to propagate from 
their sources undeflected. FuU sky coverage wifl provide maps of cosmic ray events with several 
hundred events where only a handful have ever been detected. 

Identification of a clear source of cosmic rays at extremely high energies, combined with 
observations of ultra-high energy neutrinos and gamma rays, would open a new understanding 
of cosmic accelerators. This data could solve the 50 year puzzle of the source and character 
of ultra-high energy cosmic ray accelerators. It is worth mentioning that sky maps of cosmic 
rays over the energy range observable by SWORD would provide information of strength and 
coherence length of extragalactic magnetic fields ll52ll . 

For SWORD to succeed real-time dedispersion capabilities need to be developed. Although 
the computations themselves are not complicated, this study has found that they have to be per- 
formed at a high rate. We expect that prototyping the algorithms with field-programmable gate 
arrays and implementation in application-specific integrated circuits will provide the computa- 
tional capabilities required within the power consumption constraints of the spacecraft. 
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